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From studying the time evolution of the single electron density matrix within a 
density functional tight-binding formalism we study in a fully atomistic picture the 
electronic excitation transfer between two photosynthetic pigments in real time. This 
time-dependent quantum dynamics is based on fully atomistic structural models of 
the photosynthetic pigment. We analyze the dependence of the electronic excitation 
transfer with distance and orientation between photosynthetic pigments. We compare 
the results obtained from full quantum dynamics with analytical ones, based on a 
two level system model were the interaction between the pigments is dipolar. We 
observed that even when the distance of the photosynthetic pigment is about 30 A 
the deviation of the dipolarity is of about 15 percent. 
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I. INTRODUCTION 


Photosynthesis is a natural process that begins with the absorption of sunlight by an 
arrangements of photosynthetic pigments embedded into a proteic matrix known as light 
harvesting (LH) antenna complexes. In this process, the electronic excitation (known as 
exciton) spreads and moves within and between pigments within the complex, eventualy 
reaching the reaction center where charge separation occurs^, leading to the biochemical 
energy conversion process. The quantum efficiency of the early event, namely, the absorption 
of a photon and energy transfer, is very close to 100%. A number of studies on photosynthetic 
complexes suggests that this high efficiency can be due to suppression of environmental 
decoherence of excitons within antenna complexes^ - - even a room temperature^-. Protein 
matrices in the antenna complex appear to suppress decoherence processes and allow the 
joint excitation of many photosynthetic pigments and the consistent evolution of the exciton 
for a relatively long time similar to the electronic excitation transfer (EET) time scales, 
allowing the parallel exploration of a manifold of relaxation paths during this period^ - — . 
This evidence implies that quantum coherences should be taken into account during the 
study of energy transfer in photosynthetic antenna systems. 

Understanding the physical principle underlying highly efficient energy transfer in photo¬ 
synthesis has motivated researchers to develop several theoretical models with this aim. EET 
within photosynthetic complexes is usually described by two perturbative limits depending 
on the magnitude of electron-phonon coupling. When the electronic coupling between the 
photosynthetic pigments is small, this coupling is described in terms of inter-pigment dipole- 
dipole interaction which is the lowest order nonvanishing term in the multipole expansion. 
The two assumptions underlying this approach are that vibrational relaxation within the 
excited state is much faster than energy transfer, and that the coupling to the vibrational 
modes of the bath is stronger than the coupling between the chromophores. The latter 
approach ensures that EET occurs via an incoherent hopping mechanism and shows Marko¬ 
vian behavior—In this regime the rate of EET is generally described by Forster theory^, 
where this rate can be expressed in terms of the spectral overlap of the donor emission 
spectrum and the acceptor absorption spectrum. This approach has been implemented for 
describing the physics involved in the study of EET in the coupled B800 dimer ring of the 
purple bacteria LH2 comple x 15 ^ 6 . 
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Furthermore, in the strong electronic coupling case, the excited states of the pigments mix 
and new delocalized states emerge. On short time scales, this delocalization leads to coherent 
non-Markovian dynamics in which the excitation travels as a wave until dephasing destroys 
the coherence^’ 17 . When the distance between the chromophores is similar to the spatial 
extent of the molecules, electronic wavefunctions begin to overlap and electronic exchange 
becomes important. Dexter theory includes the exchange interaction in the description of 
the EET—. However, in the case where the electron-phonon coupling is strong, neither 
Forster nor Dexter are sufficient to describe the EET. Where EET can occur on timescales 
faster than vibrational relaxations^ Redfield or Lindblad relaxation theory are usecS^ 
in order to include explicitly the coupling between chromophores and some thermal bath 
vibrational modes. 

Theoretical approaches applied to study this phenomenon are based on parametrized 
Hamiltonians containing information about inter-chromophore couplings and the excitation 
energy of individual pigments embedded in their proteic environment. Most of these pa¬ 
rameters are obtained from experimental datcS^— . Furthermore, the description of the 
inter-chromophore coupling can be obtained from quantum chemical approaches like tran¬ 
sition density cube (TDC), a three dimensional (3D) ab initio approach for the calculation 
of the EET couplings. This method gives a more realistic representation of the transition 
density of the molecules, since it takes into account the shape of the interacting molecules^. 
Madjet et al. 2& presented a numerically more efficient method for the calculation of the 
inter-chromophore coupling, known as TrEsp (transition charges from electrostatic poten¬ 
tial method). The coupling is obtained by fitting the electrostatic potential of the transition 
density on a 3D grid using atomic partial charges. 

Despite the existence of a variety of models and theoretical methods that complement 
the experimental results, to our knowledge, there is currently no method that can describe 
the photophysics of individual photosynthetic pigments as well as the dynamical evolution 
of their coupled excitations when embedded within a proteic environment from an atomistic, 
time dependent perspective^^. The main problem lies in that due to the large size of the 
antenna complex, i.e. a few thousand atoms, this entire system can not be simulated with 
the present computational methods found in literature 2 ' 29 . 

In a previous work we calculate the excitation energies as well as the transition dipole 
moments vector for a series of important photosynthetic pigments^Mi In this paper we go a 
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step further towards the comprehension of the EET mechanism. We study the full dynam¬ 
ical evolution of the one electron density matrix within a density functional tight-binding 
(DFTB) Hamiltonian in response to laser illumination. From this simulations we obtain 
the expectation value of the dipole moment, which takes into account the entire structural 
information of the molecule. From linear response theory, considering each pigment as a 
two level system (TLS) we obtain the analytical expressions of this expectations values and 
compare them with the calculated ones. As expected, the interaction between a chlorophyll 
a (Chla) dimer deviates from the dipole-dipole interaction when the distance between them 
is less or similar to the spatial extent of the wave function of each pigment, however when 
the dimer is about 30 A the interaction deviates by a 10%. This is not a minor result given 
that the models used to simulate the couplings between the photosynthetic pigments in the 
study of the energy transfer takes into account only dipolar interactions^^. We describe 
as well deviations from the dipolar coupling model as a function of the orientation of dimer 
components. 


II. COMPUTATIONAL METHOD 

The description of the electronic structure of photosynthetic pigments was carried out 
using the self consistent density functional tight-binding (SCC-DFTB) method 33 ' 31 . This 
method has been successfully applied to the description of the electronic structure of large 
molecular systems^. It is based on a second order expansion of the Khon-Sham energy 
functional around a reference density of neutral atomic species. The DFTB+ code is com¬ 
putational implementation of the DFTB method. We have used this code to obtain the 
Hamiltonian, the overlap matrix and the ground state single electron density matrix. Our 
implementation differs from that of reference 35 in that it propagates the one electron density 
matrix instead of the single particle orbitals. We have successfully applied this model to the 
calculation of chlorophyll spectra, transition dipole moment direction and magnitude. We 
obtained excellent agreement with experimental and theoretical results22^i. For the geome¬ 
tries of the photosynthetic pigments used iik&^i we applied to their ground state density 
matrix a classical sinusoidal electric held in the dipole approximation: 

H = H 0 + E 0 sin(u;Q y t) • fi (1) 
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where the held is applied in the direction and in tune with the Q y electronic transition^. 
The evolution of the system can be calculated by integrating the Liouville-Von Newmann 
equation of motion in the non-orthogonal basis: 

I = P) 

The dipole moment is calculated as: 

N 

t*(t) = e^fc(*) r i (3) 

i 

where e is the elementary charge, r* is the cartesian coordinate of atom i and is the 
Mulliken charge given by 

N 

Qi(t) = ~ (pS + Sp) vv ( 4 ) 

were Z is effective nuclear charge associated to atom i. 

To study the interaction between the chromophores within the dipole approximation, 
we consider each photosynthetic pigment as a two level system (TLS), where the difference 
between each energy level is equal to the Q y electronic transition. In the linear regime, 
where the intensity of the electric held is small, the response of the dipole moment, in 
matrix notation, to the laser perturbation is: 

j r oo ^ 

Pa(t) = -~ dr 2 _\{[P‘a('r), fip])Ep(t — r) (5) 

n Jo 8 

where a equals to x, y or z and ([/t a (r), £ip\) is the polarizability tensor that describes 
the dipole moment response in direction a to an applied electric held in direction f3. The 
expression of this tensor for a TLS is given by the following expression: 

<LMr),/V]) = j\p Qv \ 2 sin (ru Qy ) r%y Qy (6) 

where \pq v \ 2 is the transition dipole moment, c uq v = — ^ and tq v = rq y i + Tq j + TqJz 
such that the inner product of fq y is equal to 1. 

If the applied held is in the direction of the Q y electronic transition, then the expectation 
value of the dipole moment is equal to: 

2 9 r 

Pa(t) = ftE 0 \n Qy \ 2 J drsm(uj Qy T)sm(ujQ y (t-T))r% y ( 7 ) 
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FIG. 1. Schematic representation of two Chi’s molecules separated by the distance r. The molecules 
are characterized by their dipole moment, 9 is the angle between the two dipoles moments, a 
and /3 are the angles between each dipole moment and the distance vector, r. 

The solution of eq. (J7J) for long times can be approximated by the following: 

t*(t) - yl^J 2 cos(tcu Q J tf Q!/ ( 8 ) 

Suppose that the second chromophore (A) is located at f from the first one (D) and |r| 
is larger than the spatial extent of A and D, the electric field generated by the oscillating 
dipole of D is given by the following equation 

E D{t) = 1 3 (3 ■ f) • f - Mr»(t)) (9) 

47T60 r 6 

where Hoif) is the expectation value of the dipole moment of D given by eq. (J8J). e 0 is the 
vacuum permittivity and r is the distance between the Mg atoms (see Fig. [[]). 

We assumed that an electric field given by eq. ([9]) perturbs the second molecule, substi¬ 
tuting the expression of the electric field into eq. (jSJ), the expectation value of the dipole 
moment of molecule A is: 


V a(*) ^ 4 7re ^° 2 r 3 l/ ig | 2 |^| 2sin ( tuj Qy ) t 2 (cos(f3) cos(a) - i sin(/3) sin(a)^ fg y (10) 

where |ha and \hd\ 2 are the transition dipole moment of the acceptor and donor molecules, 
respectively, fg is the direction of the Q y electronic transition of the acceptor molecule, a 
and /? are the angles between each dipole moment and the distance vector, r (Fig. [[]). 
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FIG. 2. Absolute value of the dipole moment recorded as a function of time of individual molecules 
of Chi a in response to laser illumination within the time-dependent DFT description. The external 
field is applied to one of the molecules (a) while the other (b) is stimulated by the full oscillating 
field of the first, including direct electronic couplings given by matrix elements within the DFTB 
Hamiltonian. The distance between Mg atoms is 21.0 A and the molecules are aligned along the 
Q y electronic transition. 

III. RESULTS AND DISCUSSION 

As a first step towards the understanding of the interaction nature between the photo¬ 
synthetic pigments, we begin analyzing the coupling between two Chla molecules from a 
full quantum dynamical simulation without any approximations other than those implied 
in the nature of the DFTB Hamiltonian. We studied the variation of the interaction when 
distance and orientation between the molecules change independently. In both instances we 
perform a comparison between the dynamical results and the analytical ones and we show 
the percentage of deviation from the ideal dipole-dipole interaction. 

The quantum dynamical simulation starts with the perturbation of one of the Cilia’s 
with a laser type perturbation (JT|) in tune with the electronic excitation (c oq v = 1.91 eV) 
and in the direction of the Q y transition. The intensity of the perturbation was of Eq = 0.1 
mV/A in order to remain in the linear response regime. The molecules are aligned in the 
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FIG. 3. (a) Log-log plot of the dipole moment of the acceptor molecule n a and (b) the percentual 
relative error (ay) as function of the distance between corresponding magnesium atoms. The solid 
line corresponds to the analytical function of the dipole moment using as transition dipole moment 
for Chi a 1.03 e 2 A 2 —. 

Q y direction and the system was allowed to evolve in time, then we calculated the dipole 
moment of each chromophore (equation [3]). Fig. [2] describes the EET from the Chla wich is 
illuminated (donor) to the second molecule (acceptor). Fig. [ 2 ] (a) shows the time dependent 
variation of the dipole moment of the monomer at which the external held is applied. It 
can be observed that the dipole moment oscillation grows linearly with time, as can be 
expected for an applied held within the linear response regime and in the absence of any 
dissipative mechanism. Fig. [ 2 ] (b) is the time evolution of the dipole moment of the second 
molecule in response to the stimulation by the oscillating held generated by the hrst. It 
is worth noting that besides the electrostatic coupling molecules may interact via non-zero 
Hamiltonian matrix elements within the DFTB Hamiltonian. This results are in agreement 
with the analytical expressions, eq. [8] and [TO] (with a and f3 zero) derived in the previous 
section. As can be seen the oscillation of the dipole moment of the non-directly illuminated 
monomer grows quadratically in time, and depends on the distance to the hrst molecule as 


The oscillation amplitudes of the dipole moment corresponding to the acceptor molecule 









(where no laser perturbation is applied) versus the distance between the magnesium atoms, 
r, are plotted in Fig. [3] (a). The solid line is the plot of the analytical function of the dipole 
moment given by eq. ([XU . Where we have used as value of the transition dipole moment of 
the molecule D and A the value obtained for the monomeric photosynthetic pigment, this is 
1.03 e 2 A 2 —. As can be noted from the plot, for long disntances the EET can be described 
through the dipolar mechanism. However, when the distance between the molecules is 
around 40 A, the transfer of the excitation starts to deviate from the dipolar mechanism, 
and a more detailed description of the electronic coupling between pairs of pigments should 
be considered, where it must be taken into account higher order terms in the multipole 
expansion and the description of the overlap between molecular orbitals. 

Furthermore, it three regimes can be distinguished for the deviation from the dipolar 
mechanism: at long distances (r > 60 A), the slope of the line corresponding to the analytical 
function is lower than the line product of the oscillations amplitudes of the dipole moment 
calculated by the electron dynamics, this observation indicates an error by excess and can 
be attributed to a systematic error caused by the fitting of the oscillations amplitudes 
of the dipole moment with a quadratic function. This error increases with the distance 
between the photosynthetic pigments, since the magnitude of the dipole of the acceptor 
molecule decreases with r -3 . The second regime is observed at intermediate distances (17 A 
< r < 60A), where the error is by excess and is due to the absence of higher order terms in 
the multipolar expansion that describes the interaction between the molecules. Finally an 
abrupt change in the regime can be observed, where the error is by defect and is attributed 
to no taking into account the overlap between the electronic states of the two molecules in 
the multipolar expansion. This overlap allows direct charge transfer between chromophores, 
disabling EET due to the decrease of the excitation cross-section of the chromophores (\fiQ y |), 
which implies a decrease in the intensity of the oscillating electric field (E^(f)) around the 
acceptor molecule. 

This is not a minor result since the computational methods used to model the electronic 
coupling between the photosynthetic pigments in order to study the the energy transfer takes 
into account only the dipolar termr^A In photosynthetic organisms the is a remarkable 
variety of antennas, which show no apparent relation to each other in terms of structure or 
even types of pigments utilized^. We constructed a histogram of the distance between the 
Mg atoms of all dimers of photosynthetic pigments present in antenna complexes found in 
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FIG. 4. Histogram of the distance between photosynthetic pigments in the complex antennas most 
common in nature^ — 

nature for which the crystal structure has been described^'— (Fig. [4j). It can be seen that 
most of the pigments are at 30 A from each other. Therefore, if the dipole-dipole interaction 
is used to describe the interaction between the photosynthetic pigments, the results will 
show an error of around 15 % for these dimers (Fig. [31(b)) . This result shows that the use 
the dipole-dipole interaction approximation to describe their coupling is not appropriate 
at relevant interpigment distances, due both to the the influence of non-dipolar terms and 
direct electronic coupling between pigments. 

In a second step we study the dependence of the electronic excitation transfer with the 
orientation of the molecules (equation HO]). For this purpose a continuous laser was applied 
to one Chla in tune with the Q y electronic excitation and the dipole moment of the acceptor 
Chla located at 20 A was computed. Several simulations were performed by varying the 
orientation of the acceptor molecule, and Fig. [5] shows the orientation dependence for the 
dipole moment of thhis molecule. It can be seen that this variation has the shape of a 
cosine function. When the two pigments are oriented perpendicular to each other (/3 = 7t/2 
rad) the dipole moment is zero, indicating that the second molecule is not excited and the 
maximum excitation is reached when the two pigments are aligned with the Q y electronic 
transition (/3 = 0 rad and (3 = tt rad). 
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FIG. 5. Dependence of the dipole moment of the acceptor molecule, with the orientation of 
this dipole moment with respect /3. the distance between the two Mg atoms of the chromophores 
is 20 A. 

To analyze in more detail the deviation from dipolar behavior, Fig. |6] shows the depen¬ 
dence of the dipole moment of the acceptor molecule with the | cos(/3)| for different distances 
between the Mg atoms (upper panel) and the percentage relative error between the analyt¬ 
ical results (eq. (ITOj) ) and the results obtained from eq. (J3J) (lower panel). If the interaction 
is dipolar the dipole moment of the acceptor molecule varies linearly with respect | cos(/3)|. 
It can be observed that there is a greater deviation from the ideal behavior for r = 15 A 
and as the distance between the molecules starts to increase the deviations decreases. From 
the plot of the relative error we can observed that the behavior is not equal when the angles 
are supplementary, this is due to the fact that the molecules are not symmetric. 

When the distance between the pigments is 15 or 20 A the largest percentage of de¬ 
viation from dipolarity is when the pigments are aligned and the lower deviation is when 
the acceptor molecule is oriented at 7t/4 rad, this is because this orientation allows closer 
approximation between the pigments. When the distance between the molecules is 30 and 
50 A the deviation from the dipolarity begins to decrease and when they are aligned the 
deviation is of 5%. Once again, it is worth noting that this is not a minor result, since most 
of the photosynthetic pigments found in nature (Fig. [3J are located 30 A from each other, 
in some antenna complexes the distances are even lower. Also, most of the models used to 
study the EET in these systems are based in the exciton Hamiltonian, where the electronic 
coupling between the pigments is essentially dipolar—*—therefore, if the dipolar interac- 
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FIG. 6. Dependence of the dipole moment of the acceptor molecule (upper panel) and the per¬ 
centage of the relative error (lower panel) with | cos(/3)|. 

tion is used to describe the interaction between the chromophores the associated error will 
of the order of 15% for chromophores at distances between 20 and 30A. Finally, it can be 
noted that when the orientation of the acceptor molecule is about 7 t/ 6 rad the deviation 
is bigger than when the two chromophores are aligned. This deviation from the dipolar 
interaction is due to systematic errors made when the oscillations amplitudes of the dipole 
moment of the acceptor molecule is fitted to a quadratic function. The error is greater in 
this orientation because the electronic coupling and the magnitude of the dipole moment of 
the acceptor molecule are lower than those obtained for the collinear orientation. 


IV. CONCLUSIONS 

In this work we describe the energy transfer between two Chla molecules in real time on 
the basis of a fully quantum time dependent model as a function of of their distance and 
orientation. Results for the time dependence of the dipole moment of the two components of 
the dimer are obtained when one of the pigments is illuminated with a time dependent held in 
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tune with its Q y excitation. Results for the time dependent variation of the dipole moments 
of the illuminated and non-illuminated monomer are interpreted in terms of a simple two 
level system scheme. The growth of the dipole signal with time that results from the full 
dynamical evolution of the dimer density matrix within the DFTB model is the compared 
to that expected from the consideration of a dipole-dipole coupling between pigments. Our 
results indicate that significant differences exist between the simple dipole description an 
the DFTB Hamiltonian which are more important for dimers placed at distances of about 
30 A or closer, which are not uncommon in antenna complexes found in nature. 
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